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ABSTRACT 

Based on our long (~ 300 ks) 2007 XMM-Newton observation of the Seyfcrt galaxy NGC 1365, 
we report here on the spectral and timing behaviour of two ultraluminous X-ray sources, which had 
previously reached isotropic X-ray luminosities Lx ~ 4 x 10 40 erg s" 1 (0.3-10 keV band). In 2007, 
they were in a lower state (Lx ~ 5 x 10 39 erg s _1 , and Lx ~ 1.5 x 10 39 erg s _1 for XI and X2, 
respectively). Their X-ray spectra were dominated by power-laws with photon indices r ra 1.8 and 
r ra 1.2, respectively. Thus, their spectra were similar to those at their outburst peaks. Both sources 
have been seen to vary by a factor of 20 in luminosity over the years, but their spectra are always 
dominated by a hard power-law; unlike most stellar-mass BHs, they have never been found in a 
canonical high/soft state dominated by a standard disk. The lack of a canonical high/soft state 
seems to be a common feature of ULXs. We speculate that the different kind of donor star and/or 
a persistently super-Eddington accretion rate during their outbursts may prevent accretion flows in 
ULXs from settling into steady standard disks. 

Subject headings: X-rays: binaries — X-rays: individual (NGC 1365) — black hole physics 



1. INTRODUCTION 

One of the hallmarks of accretion onto stellar-mass 
black holes (BHs) is a pattern of state transitions, ini- 
tially discovered from X-ray spectral studies (starting 
from Cygnus X-l: Tananbaum et al. (1972)) and later 
understood as fundamental changes in the structure 
of the accretion flow, and recognized also in the ra- 
dio and optical behavior. The "canonical" scheme of 
BH accretion (Remillard & McClintock 2006; Fender, 
Belloni & Gallo 2004; Esin et al. 1998) consists of at 
least three main states: low/hard, high/soft and very 
high. The low/hard state is dominated by a power- 
law spectrum in the X-ray band, probably produced by 
inverse- Compton scattering of soft seed photons (Sun- 
yaev & Titarchuk 1980). The location and geometry of 
the scattering region is still unclear: it has been sug- 
gested that it could be a hot corona above the cold 
disk (Haardt & Maraschi 1993; Poutanen & Svensson 
1996), an advection-dominated accretion flow (ADAF) 
(Narayan & Yi 1994; Esin, McClintock k, Narayan 1997), 
a fast jet or the base-of-the-jet (Fender, Belloni & 
Gallo 2004), or a hot, centrifugally-supported boundary 
layer or post-shock region due to sub-Keplerian inflows 
(Chakrabarti & Titarchuk 1995; Chakrabarti 1997). The 
high/soft state is dominated by thermal X-ray emission; 
most of the accretion power is radiated by an optically- 
thick, geometrically-thin disk (Shakura & Sunyaev 1973). 
The very high state, observed when a source approaches 
its Eddington luminosity, is still poorly understood; it is 
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again dominated by a power-law component in the X- 
ray band (usually steeper than in the low/hard state, 
Remillard & McClintock (2006)), together with a ther- 
mal disk component that becomes relatively less impor- 
tant as the luminosity increases. The power-law com- 
ponent is thought to be produced by the Comptoniza- 
tion of soft X-ray photons from the inner disk region 
(Done & Kubota 2006). The spectrum of the disk com- 
ponent itself may be heavily modified by electron scat- 
tering, energy advection and radiation trapping, in the 
very high state (slim disk model: Watarai, Mizuno & 
Mineshige (2001)). In the simplest approximation, an 
accreting BH is in the radiatively-inefncient low/hard 
state when the mass accretion rate is less than a few per 
cent of the Eddington accretion rate (i.e., the accretion 

parameter m = M/M Edd ra (o.lc 2 M) /L Edd < 0.01). 

At higher accretion rates (0.01 < m < 0.5), accreting 
BHs are expected to switch to the radiatively-efficient, 
disk-dominated high/soft state and then to the very 
high state, until the standard disk becomes unstable (at 
ra > 1). Another fundamental property, revealed by ra- 
dio observations, is the presence in the low/hard state 
of a steady jet which may carry most of the outgoing 
accretion power; the jet is suppressed in the high/soft 
state (Merloni, Heinz & di Matteo 2003; Fender, Belloni 
& Gallo 2004) ; strong X-ray and radio flares are instead 
observed in the very high state (Remillard & McClintock 
2006). 

One of the interesting properties of the standard disk 
spectrum in the high/soft state is that we expect a sim- 
ple relation between the disk luminosity, the disk peak 
temperature, and the BH mass (Makishima et al. 1996, 
2000): i disk ra AirR^aT^ ~ Af| H 7£. In this formula, 
the inner-disk radius i?; n ra 6i? g for a non-rotating BH, 
where the gravitational radius R g = Mbh in geometric 
units. This relation has been observationally confirmed 
for many Galactic BHs, and the BH masses inferred with 
this method are generally in agreement (within a factor of 
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2) with direct kinematic measurements. Another useful 
property of the high/soft state is the simple parameter- 
ization of disk luminosity and temperature with the ac- 
cretion rate: for a fixed BH mass, Ldisk ~ ih t T- m ~ rh 1 ^ 4 
(Shakura & Sunyaev 1973). 

The canonical state classification has also been ap- 
plied to AGN, justified by the self-similar nature of 
many physical processes near a BH horizon, provided 
that characteristic sizes, timescales and disk tempera- 
tures are properly rescaled by the BH mass (Maccarone, 
Gallo & Fender 2003; Merloni, Heinz & di Matteo 2003; 
Jester 2005; Markowitz & Uttley 2005). In this scheme, 
low-luminosity AGN and LINERs may represent the 
low/hard state, typical AGN and Seyferts may be in the 
high/soft state, and the most powerful quasars and FR-II 
radio galaxies may be in the very high state. 

Ultraluminous X-ray sources (ULXs) are accreting 
BHs with observed fluxes corresponding to isotropic X- 
ray luminosities up to ~ a few 10 40 erg s _1 , well above 
the Eddington limit for the kind of stellar BHs known in 
our Local Group (Roberts (2007) for a recent review). It 
is still widely debated whether this is due to extremely 
heavy stellar BHs, with masses ~ 30-70M Q (Pakull & 
Mirioni 2002), or intermediate-mass BHs (Miller & Col- 
bert 2004), or collimation of the outgoing radiation (King 
et al. 2001; King 2008), or super-Eddington luminosity 
(Begelman 2002, 2006; Ohsuga & Mineshige 2007); or 
perhaps a combination of some of those factors. In other 
words, it is still unclear whether ULXs represent a differ- 
ent physical class of accreting BHs (perhaps with their 
own sequence of canonical states), or a different accretion 
state of "ordinary" stellar-mass BHs. 

In the absence of direct kinematic measurements of BH 
masses in ULXs (mainly because of the faintness of their 
optical counterparts, at typical distances > a few Mpc), 
X-ray spectroscopy could in principle provide strong con- 
straints, by analogy with the properties of canonical spec- 
tral states in stellar-mass BHs. For example, if we ob- 
serve a ULX in different spectral states over time, we 
could estimate its luminosity at the transition between 
low/hard and high/soft state, or the maximum luminos- 
ity in the low/hard state, and the fitted disk luminos- 
ity and temperature in the high/soft state. Such argu- 
ments were used for example by Winter, Mushotzky & 
Reynolds (2006) and Miller & Colbert (2004), in support 
of BHs with masses > 10 3 M Q . However, this method 
relies on the assumption that ULXs follow the same pat- 
tern of canonical spectral states; specifically, that they 
will switch from the low/hard to the high/soft state at a 
luminosity ~ a few per cent of LecM, an( i that they are 
dominated by a standard disk in the high/soft state. 

Until now, it has been difficult to test these assump- 
tions. The main reason why we still know very little 
about spectral states transitions in ULXs is that none of 
these sources has had the kind of daily or weekly X-ray 
monitoring that is possible for Galactic BHs. Even the 
brightest ULXs are usually observed only sporadically, 
every few months or years. This makes it difficult for us 
to understand the sequence and timescale of their spec- 
tral state transitions. Nonetheless, for a few ULXs there 
are now enough datapoints over the years, so that we can 
try and compare their behaviour with those of transient 
Galactic BHs. 



Fig. 1. — EPIC-M.OS true-color image of the central region of 
NGC 1365. Red = 0.3-1 keV; green = 1-2 keV; blue = 2-10 keV. 
The two ULXs discussed in this paper are marked as "1" and "2" . 

NGC 1365 XI is one of such rare sources, detected by 
ROSAT, ASCA, Chandra, Swift and XMM-Newton at 
various levels of activity between 1993 and 2007 (Ko- 
mossa & Schulz 1998; Soria et al. 2007; Strateva & Ko- 
mossa 2008), with peak luminosities w 3-4 x 10 40 erg 
s" 1 in 1995 and 2006 (Soria et al. 2007). The nearby 
transient ULX X2 also reached a luminosity «4x 10 40 
erg s _1 in 2006 (Strateva & Komossa 2008), but was 
much fainter or undetected during the earlier observa- 
tions. Their host galaxy is the most luminous spiral in 
the Fornax cluster, [type SBb(s)I: Sandage & Tammann 
(1981)] with a total gas mass ss 3 x 10 10 Mq, a kinematic 
mass « 3.6 x 10 n Mo, and a star formation rate w 10M© 
yr" 1 (Lindblad 1999; Roussel et al. 2001). The Cepheid 
distance to NGC 1365 is 19 ± 1 Mpc (Ferrarese et al. 
2000); at that distance, 1 arcsec = 92 pc. In this paper, 
we report the results of our ps 300-ks XMM-Newton ob- 
servation, split over three orbits between 2007 June 30 
and 2007 July 05 (Table 1 and Section 2 for details). We 
interpret the X-ray properties of the two ULXs (in par- 
ticular, their hard power-law spectra) in the context of 
their long-term behaviour. 

2. OBSERVATIONS AND DATA ANALYSIS 

XMM-Newton observed NGC 1365 on 2007 June 30- 
July 05 (Table 1); the Principal Investigator was Guido 
Risaliti. The prime instrument was the European Pho- 
ton Imaging Camera (EPIC), in the Full Window mode, 
Medium filter. See Figure 1 for a true-colour image of 
the combined MOS dataset. The primary target was the 
highly variable Seyfert nucleus, on which we report else- 
where (Risaliti et al. 2008, in prep.). We used the XMM- 
Newton Science Analysis System (SAS) version 7.1.0 to 
process and filter the event files and extract lightcurves 
and spectra. Several background flares occurred during 
the three observations. We tried various rejection thresh- 
olds to increase the signal-to-noise for the non-nuclear 
sources; in particular, to improve the signal-to-noise of 
XI and X2 at photon energies > 5 keV (the two sources 
have similar fluxes in that energy range). Instead, the 
nucleus is brighter, and therefore much less affected by 
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TABLE 1 

Log of the 2007 XMM-Newton observations of NGC 1365 



Revolution Obs ID Instrument Start time End time Live time (ks) Clean GTI (ks) 



1384 0505140201 MOS 2007-06-30 at 07:08:45 2007-07-01 at 18:51:55 126.7 96.6 

pn 2007-06-30 at 07:31:03 2007-07-01 at 18:48:52 110.8 78.4 

1385 0505140401 MOS 2007-07-02 at 07:07:24 2007-07-03 at 18:23:52 124.7 123.3 

pn 2007-07-02 at 07:39:43 2007-07-03 at 18:35:02 109.2 107.9 

1386 0505140501 MOS 2007-07-04 at 08:08:58 2007-07-05 at 18:27:14 121.7 92.8 

pn 2007-07-04 at 08:36:56 2007-07-05 at 18:24:16 105.8 55.9 



the flares, except for a particularly strong, short one at 
the end of the second observation. In the end, for our 
spectral analysis we chose to retain good-time-intervals 
of 78.4 ks, 107.9 ks, and 55.9 ks for the EPIC-pn in the 
three observations, and 96.6 ks, 123.3 ks, and 92.8 ks 
for the EPIC-MOS (Table 1). For the timing analysis, 
we used the whole live exposure time, although of course 
the intervals with high or flaring background give lower 
signal-to-noise bins in the lightcurve. 

The extraction of spectra and lightcurves for XI was 
made difficult by the unfortunate location of the source 
on a chip gap in the pn, leading to a loss of about half the 
counts in every exposure. The ancillary response func- 
tion is supposed to take into account the reduced effec- 
tive area of the source region; however, because the point 
spread functions of hard and soft photons are slightly dif- 
ferent, a chip gap may in principle introduce a change in 
the spectral shape. Fortunately, XI is not on a chip gap 
in the MOS images. We fitted separately and compared 
the MOS and pn spectra, and verified that they are con- 
sistent with each other. Thus, we conclude that the chip 
gap may have decreased the signal-to-noise achievable in 
principle, but has not significantly affected the spectral 
results. 

The radius of the source extraction regions was 20" 
for XI (the brightest non-nuclear source), and only 8" 
for X2, to reduce contamination from a nearby source of 
similar brightness (X2 was much fainter than in 2006). 
Background extraction regions were chosen around the 
source regions, in a suitable way to avoid contamina- 
tion. After building response and ancillary response files 
with rmfgen and arfgen we used XSPEC (Arnaud 1996) for 
spectral analysis. To improve the signal-to-noise ratio, 
we co-added the EPIC pn and MOS spectra from each 
observation, and for all the observations, with suitably 
averaged response functions, using the method of Page 
et al. (2003). 

3. MAIN RESULTS 

3.1. Time variability of XI 

We extracted background-subtracted MOS and pn 
lightcurves of XI in the three observations, and analyzed 
them with standard FTOOLS tasks, with various choices 
of binning. We did not find any significant variability 
between the three observations, nor within a single ob- 
servation. For example, for the 2007 July 02-03 EPIC-pn 
observation (Rev. 1385), which has the highest signal-to- 
noise, we find a \ 2 probability of constancy > 0.9 (Figure 
2). We do not find any significant features in the power 
density spectra, either. 




2x10" 4x10 4 6x10" 8x10" 10 5 1.2x10 5 

Time (s) 

Fig. 2.— EPIC-pn lightcurve of XI (binned to 1000 s) from the 
second XMM-Newton observation of 2007 (Rev. 1385). 



Over the past fourteen years, XI has shown strong 
variability, with two major outbursts detected by ASCA 
in 1995 and Chandra in 2006 (Figure 3). In addi- 
tion to the luminosity datapoints already summarized 
in Soria et al. (2007), we have also analyzed a series 
of short Swift monitoring observations, carried out on 
2006 July 21-24 (three months after the Chandra obser- 
vation), for a total exposure time of 19.4 ks. We de- 
tect 23 net counts with the X-Ray Telescope (count rate 
w 0.0012 count s" 1 ), corresponding to an emitted lumi- 
nosity w (2.3±0.5) x 10 39 erg s" 1 in the 0.3-10 keV band 
(Figure 3) , assuming a power-law model with photon in- 
dex 1.8 and iV H = 5 x 10 20 cm" 2 . 

3.2. Spectral analysis: XI is a power-law ULX 

Firstly, we fitted and compared the MOS spectra of 
XI from each XMM-Newton observation during the three 
consecutive revolutions of 2007 June 30 - July 05. We 
find that the three spectra are consistent with each other 
with regard to observed flux, absorbing column density 
and spectral shape (Figure 4) . All three spectra are well 
fitted by an absorbed power-law model with total ab- 
sorbing column density ~ 5 x 10 20 cm -2 and power-law 
photon indices « 1.7-1.8, consistent within the 90% con- 
fidence limit of each dataset. Hence, we are justified in 
adding the EPIC spectral data from all three observa- 
tions to increase the signal-to-noise ratio. 

The combined XMM-Newton /EPIC spectrum is well 
fitted (xl = 1.03(156.2/152)) by a simple power-law 
model (see the top panel of Figure 5, and Table 2 for 
the best-fitting parameters). From this fit, we estimate 
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Fig. 3. — Long-term X-ray lightcurve of XI: unabsorbed lumi- 
nosities in the 0.3—10 keV band. Green filled circles are the lu- 
minosities extrapolated from ROSAT; open magenta circles from 
ASCA; red squares from Chandra; blue triangles from XMM- 
Newton; an open cyan circle from Swift. For details on the ob- 
servations from 1993 to early 2006, see Soria et al. (2007). The 
peak ASCA luminosity PS 3-4 X 10 40 erg s _1 is based on our own 
re-analysis of the data (Soria et al. 2007), and is about half of the 
value often quoted in the literature. 



an unabsorbed, isotropic luminosity Lx ~ 5.3 x 10 39 erg 
s~ 4 in the 0.3-10 keV band. Some ULXs often show 
a curvature in their X-ray spectra (see Section 4), to 
the point that it is sometimes difficult to distinguish 
between a power-law, a cut-off power-law and a disk- 
blackbody model, ft is clearly not the case here: we 
show the best-fitting disk-blackbody model in the bot- 
tom panel of Figure 5 {xl > 2), as a graphic evidence 
of the lack of spectral curvature. For the same rea- 
son, even a bremsstrahlung model, although slightly less 
curved than a disk-blackbody, does not fit the data well; 
the best-fitting model has xl = 1-13(172.3/152), for 
kT = (5.6 ± 0.6) keV. 

A power-law is the best-fitting 1-component phe- 
nomenological model, and it makes no assumptions on 
the underlying physical process. It is of course possible 
to fit the same spectrum equally well with other, more 
complex physical models, provided that they look exactly 
like a power-law in the observed band. For example, we 
tried the Comptonization model comptt and found that 
we could get equally good fits by fixing the seed-photon 
temperature at any value fcTb < 0.1 keV, if the coronal 
temperature kT e > 2 keV and the optical depth t > 0.5. 
The best-fitting comptt model has xl = 1-03(155.2/151) 
for kT e w 2.9 keV and r « 6.2 (warm, optically-thick 
corona). But we also get a statistically equivalent fit 
{xl = 1-04(157.8/152)) for example for kT e > 50 keV 
and r = 0.5 (hot, optically-thin corona). Both models 
are indistinguishable from a power-law, with the avail- 
able data. 

We then tried adding a diskbb component (soft ex- 
cess) to the power-law: the best-fitting model has xl — 
1.02(153.2/150) for kT in = 0.35±g;^ keV and photon 
index T = 1.66^q'^, but this is not a statistically- 
significant improvement with respect to the simple 
power-law model. From the upper limit to the diskbb 
normalization, we find that, if there is an additional disk- 
blackbody component, it does not contribute more than 



« 3.6 x 10^ erg s" 1 in the 0.3-10 keV band (« 5.4 x 10 
erg s _1 bolometric), that is < 7% of the power-law con- 
tribution. 

Finally, we looked for a possible break or cut-off in 
the power-law at high energies. An exponential cut- 
off does not improve the fit, because it introduces ex- 
cessive curvature in the model. A break at the energy 
E h = 7.5 ± 0.9 keV gives xl = 1-00(150.3/150), which 
is, formally, a marginally significant improvement from a 
simple power-law fit, and confirms the visual impression 
of a spectral downturn (Figure 5, top panel). However, 
given the low signal-to-noise level at those channel ener- 
gies, and the uncertainties in the background subtraction 
(contaminated by the strong AGN emission above rs 6 
keV), we cannot attribute any physical meaning to this 
possible break. Nonetheless, we recall that there was also 
marginal evidence of a break at sa 6 keV in one of the 
XMM-Newton observations of 2004 (Soria et al. 2007). 

The power- law photon index T k 1.8 (Table 2) is 
consistent with the indices fitted to the Chandra spec- 
tra during the 2006 April outburst (see Table 2 for the 
best-fitting parameters of the combined 2006 April 12 15 
spectrum), and also with the indices fitted to the ASCA, 
Chandra and XMM-Newton spectra over 1995-2003. In 
the 2004 XMM-Newton spectra, instead, XI showed a 
soft excess (inner disk temperature kT- ln w 0.3-0.4 keV) 
contributing about 25-30 per cent of the X-ray luminos- 
ity, in addition to a dominant (and harder) power-law 
component with a high-energy break. At that time, the 
unabsorbed luminosity was rs 10 40 erg s" 1 (Soria et al. 
2007). 

The column density in the 2007 spectrum is consistent 
with the values previously fitted to the XMM-Newton 
spectra when XI had X-ray luminosities < 10 40 erg 
s _1 . The neutral column density was instead a few 
times higher during the 2006 outburst, when a single- 
component power-law model was used (Soria et al. 2007). 
As a further check that the decrease in column den- 
sity between 2006 and 2007 is real and not due to the 
appearance of an underlying soft component, we fitted 
the 2007 spectrum with an additional diskbb component 
(leaving its temperature and normalization as free pa- 
rameters) and Ah fixed at the 2006 value. We obtain 
xl = 1.15(173.9/151) which is significantly worse than 
the low-column best fit. Besides, there was also strong 
ionized absorption at the peak of the 2006 outburst (So- 
ria et al. 2007), which is absent from the 2007 spectra. 

We also re-examined the 2006 data to test the opposite 
scenario: the possibility that the column density in 2006 
was as low as in 2007 or at any other epoch. We started 
from the best-fitting model {xl = 0.72(105.5/146)) to 
the 2006 outburst summarized in Table 7 of Soria et 
al. (2007), but this time we imposed A H = 4.5 x 10 20 
cm~ 2 . A single power-law model is no longer acceptable 
{xl = 1.06(155.5/147)), because it does not reproduce 
the low-energy downturn. A 1ow-Ah disk-blackbody 
model is statistically worse than the best-fitting model 
{xl = 0.76(113.2/146)): it can reproduce the low-energy 
downturn but, in that case, it overpredicts high-energy 
curvature; the best-fitting kT- ln m 1.4 keV (a characteris- 
tic temperature that makes it difficult for Chandra/ ACIS 
to distinguish between an absorbed power-law and a disk- 
blackbody model, because of its limited spectral cover- 
age). On the other hand, low- Ah Comptonization mod- 
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TABLE 2 

Best-fitting parameters for the combined 2007 
XMM-Newton/EPIC spectrum of XI, compared with the 
2006 peak outburst spectrum from Chandra data. Spectral 
MODEL: phabs*phabs*power- law. ERRORS ARE 90% CONFIDENCE 
LEVELS FOR 1 INTERESTING PARAMETER (A% 2 = 2.7). 



TABLE 3 

Best-fitting parameters for the combined 2007 
XMM-Newton/EPIC spectrum of X2. compared with the 

2006 PEAK OUTBURST SPECTRUM FROM Chandra DATA. SPECTRAL 
MODEL: phabs*phabs*power- law. ERRORS ARE 90% CONFIDENCE 
LEVELS FOR 1 INTERESTING PARAMETER (Ax 2 = 2.7). 



Parameter 


zuuo Lstianara, value 


zuu/ AiWiw-iv ewton value 


AfH.Gal" 

r b 


1.3 x 10 20 
l.5±g;| x 10 21 

1 7 4+0.i2 
L - '^-0.11 

9.8±\\l x 10- 5 


1.3 x 10 20 
4.3±\i x 10 20 
i sn+ 04 

1 - 8U -0.05 

1 - 9 -to!i x 10-5 


x 2 /dof 


0.93(85.4/92) 


1.03(156.2/152) 


/o.3-10 d 
£o.3-10 6 


5.2±° l x 10- 13 
2.8± 9 ; 3 x 10 40 


l.l±°-\ x 10" 13 
5.3± ,; 3 x 10 39 



a From Kalberla et al. (2005). Units of cm" 2 . 
b Photon index. 

c Units of photons kcV~ 1 cm -2 s — 1 , at 1 keV. 

d Observed flux in the 0.3—10 keV band; units of erg cm -2 s — 1 . 

°Unabsorbcd luminosity in the 0.3-10 keV band; units of erg s -1 . 

els provide statistically-equivalent fits to the high-A^H 
power-law model, if the seed photon temperature w 0.3 
keV. For example, using the Comptonization model bmc, 
we obtain xl = 0.72(104.5/145) with the best-fitting 
seed photon temperature fcTo = 0.27+o'o4 keV. The same 
Comptonization models applied to the 2007 spectra re- 
quire a lower seed photon temperature, < 0.1 keV. As 
for the ionized absorption component seen at the peak of 
the 2006 outburst (Soria et al. 2007), we have not found 
it possible to eliminate it or replace it with any other 
simple spectral component. 

3.3. X2: a transient, hard power-law ULX 

The transient ULX X2 (Figure 1) was serendipitously 
discovered during the 2006 Chandra observations (Strat- 
eva & Komossa 2008) 6 . It underwent an outburst with 
a rapid rise time (a factor of 7 in two days; Strateva & 
Komossa (2008)) at the same time as the unrelated ULX 
XI was also in outburst. In our 2007 XMM-Newton ob- 
servations, the source was much fainter, and much less 
variable; the observed flux is the same in the three obser- 
vations, within lcr. Remarkably, the slope of the X-ray 
spectrum was similar in 2006 and 2007, well fitted by a 
hard power-law with L « 1.1—1.2 (Table 3 and Figure 6). 
The column density in 2007 was « 10 21 cm -2 , an order 
of magnitude less than at the peak of the 2006 outburst. 
The unabsorbed isotropic luminosity reached w 3.7 x 10 40 
erg s" 1 in 2006, and was down to ~ 1.5 x 10 39 erg s _1 in 
2007. 

As for XI, we tried fitting the 2007 XMM-Newton 
spectrum of X2 with alternative models, even though 
the lower signal-to-noise makes multi-component mod- 
els even less constrained. A diskbb model provides a 

6 It is not the same source designated as X2 in Fig. 2 of Soria et 
al. (2007). 



Parameter 2006 Chandra Value 2007 XMM-Newton Value 


Ari,Gai a 


1.3 x 10 20 


1.3 x 10 2() 




8.1± 2 j x 10 21 


1.2±g;J x 10 21 


r b 


1 oq+0.25 
1 - ei> -0.19 


, 1O+0.09 
1 - L - i -Q.W 




7.6l 3 ; 2 x IO" 5 


2.6±° Q 6 5 x 10- 6 


X 2 /dof 


0.70(23.2/33) 


1.12(110.8/99) 


/o.3-10 d 


6.8+3 2 x 10_13 


3.2+. ,; 4 x 10- 14 


io.3-10 6 


3.7t°; 3 x 10 40 


l.5±g;i x 10 39 


a From Kalberla et al. (2005). Units of cm~ 2 . 
b Photon index. 



c Units of photons keV -1 cm -2 s — 1 , at 1 keV. 

d Observed flux in the 0.3-10 keV band; units of erg cm -2 s — 1 . 

°Unabsorbed luminosity in the 0.3-10 keV band; units of erg s -1 . 

formally acceptable fit (x 2 = 1.15(113.6/99), compared 
with xl = 1.12(110.8/99) for the power-law model). 
But the best-fitting temperature kT ln = 3.48iJ;|g keV 
is unphysically high for a standard disk, and the nor- 
malization corresponds to an unphysically small inner 
radius r; n (cost?) ' 5 7 km. A bremsstrahlung model 
provides a formally acceptable fit (x 2 = 1.14(112.6/99)) 
for kT > 48 keV. The Comptonization model comptt 
with seed photon temperature fixed at fcTo = 0.1 keV 
and optically thin corona (r = 0.5) gives an equivalent 
fit {xl = 1.12(111.3/99)) for kT e > 250 keV. For an 
optically-thick corona (r = 10), the best-fitting model 
has kT e = 2.73±g"Jl keV (x 2 = 1.11(109.7/99)). An 
exponential cutoff does not improve the power-law fit 
(x 2 = 1.12(110.8/99)), and the cutoff is constrained to 
be at energies > 5.8 keV. A break does not improve the 
power-law fit, either. Adding a soft disk-blackbody com- 
ponent to the power-law does not improve the fit: the 
best fitting model has a null normalization for diskbb, 
and an upper limit of w 25% for the disk-blackbody 
contribution. In short, as was the case for XI, also 
X2 can be fitted by various physical models, as long as 
they look exactly like a hard power law in the observed 
band. We then fixed the absorbing column density to 
the 2006 value (Table 3) and verified that it is not con- 
sistent with the 2007 data, even with the addition of 
an unconstrained disk-blackbody component to compen- 
sate for that. The best-fitting high-absorption model has 
X 2 = 1.25(122.5/98). Finally, we re-examined the 2006 
data, and tried fitting them with the same low absorption 
as in 2007. As for the case of XI, we find that Comp- 
tonization models with low absorption but high seed- 
photon blackbody temperature (fcT w 0.9 keV) give 
fits that are only slightly worse than the high-absorption 
models. However, such scenario is too contrived, and un- 
seen in any other accreting source. Hence, we conclude 
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Fig. 4. — Individual EP/C-MOS spectra of XI during the three 
observations in 2007, and respective \ 2 residuals, for absorbed 
power-law fits. Black datatpoints and residuals are for Rev. 1384; 
red ones for Rev. 1385; green ones for Rev. 1386. The three spectra 
are statistically identical in flux normalization and spectral slope. 

that the decrease in absorption at lower luminosities is 
physically significant for X2. 

In previous years, X2 was not seen by ASCA. It was 
claimed (Strateva & Komossa 2008) that it may have 
been marginally detected by ROSAT/RRl in 1995, at 
a count rate corresponding to a 0.3-10 keV luminosity 
w 4 x 10 39 erg s _1 , if we assume the spectral parame- 
ters of the 2007 XMM-Newton spectrum; however, this 
detection is doubtful, because of likely source confusion. 
X2 was also in a low state (Lx < 3 x 10 38 erg s _1 ) dur- 
ing the 2002 Chandra observation. It may be marginally 
detected (but confused with nearby, brighter sources) in 
the 2004 XMM-Newton observations, at a luminosity ~ 
a few 10 38 erg s _1 . It is not detected (Lx ^ 10 39 erg 
s" 1 ) in the Swift observations of 2006 July. 

4. DISCUSSION AND SPECULATIONS 
4.1. Power-law versus thermal disk spectra 

The strongly star-forming spiral galaxy NGC 1365 con- 
tains many X-ray sources that have exceeded the ULX 
threshold at some stage in the past few years. Here we 
have focused on two sources that showed recent outbursts 
and great changes in flux; for one of them (XI) there is a 
long sequence of X-ray observations dating back to 1993. 
The X-ray spectra of XI and X2 appear always domi- 
nated by broad, power-law-like components with photon 
indices T « 1.8 ± 0.1 and T rj 1.2 ± 0.1, respectively, 
both in their higher and lower states (even though their 
luminosities varied by a factor of 20). For XI, some- 
times (XMM-Newton observations from 2004, Soria et al. 
(2007)) there is also a cool, thermal soft component in 
addition to the dominant, hard power-law component. 
For both sources, a significant difference in the 0.3-10 
keV band bewteen spectra taken at higher and lower lu- 
minosity is the higher intrinsic column density in the 
higher state. We do not have enough elements to deter- 
mine whether this is directly related to a higher mass 
inflow or outflow rate, or both. 

Based on their phenomenological spectral appearance, 
one may be tempted to classify both ULXs in the power- 
law-dominated low/hard state, using the canonical def- 
inition of BH spectral states. Several other luminous 



ULXs (with L x < 3 x 10 39 erg s _1 ) are well fitted by 
a single power-law of photon index 1 < T < 2 (Berghea 
et al. 2008): such slopes are consistent with and some- 
times even harder than those measured in the low/hard 
state of Galactic BHs. However, the first hint that the 
accretion state of XI and X2 may not be the canonical 
low/hard state comes from their long-term variability: 
both sources remained in the same hard state over the 
years, despite their luminosity changes. This is different 
from the behaviour observed in most stellar-mass Galac- 
tic BHs, which recurrently switch between the low/hard 
and high/soft states. In particular, in most Galactic 
soft X-ray transients, outbursts tend to start in the hard 
state, but rapidly switch to the soft state (on a timescale 
of few days, during which the disk builds up); then, 
such sources spend the majority of their bright time (a 
few weeks or months) in the high/soft state, dominated 
by standard disk-blackbody emission with L ~ and 
0.1 < L < L Edd . 

This does not appear to be the case for ULXs: none of 
the most luminous or best studied sources has ever been 
unequivocally found in a canonical high/soft state. Some 
ULXs do have a convex spectrum (Makishima et al. 2000; 
Feng & Kaaret 2005; Stobbart et al. 2006; Tsunoda et 
al. 2007; Mizuno et al. 2007; Makishima 2007; Miyawaki 
et al. 2008), with various degrees of curvature between a 
straight power-law and a disk-blackbody. But if they are 
formally fitted with a disk-blackbody model, their peak 
color temperature (fcT; n ~ 1.5-2.5 keV) is much higher 
than expected for a standard disk, and implies a luminos- 
ity L > LEdd, inconsistent with the standard-disk model. 
One possible interpretation is that the emission is from 
a hot, super-Eddington slim disk (Watarai, Mizuno & 
Mineshige 2001; Makishima 2007). Another interpreta- 
tion is that it is Comptonized emission from a warm, 
optically-thick corona covering the inner part of the disk 
(Roberts 2007) . A second subgroup of ULXs (including 
NGC 1365 XI in 2004) sometimes show a thermal soft 
excess (kT ln rj 0.15 keV) in addition to the power-law 
component; but the soft excess contributes only < 30% 
of the luminosity (Stobbart et al. 2006). Again, this 
is not evidence of a canonical high/soft state, which is 
by definition dominated by the disk. In fact, it is not 
uncommon also for Galactic BHs in the low/hard state 
to have a non-dominant soft excess from the truncated 
outer disk in addition to the hard power-law component 
(Remillard & McClintock 2006). In summary, although 
there are two subgroups of ULXs with possible evidence 
of thermal components in their X-ray spectra, neither 
case corresponds to the canonical high/soft state. 

4.2. A different accretion state or a different class of 

BHs? 

Based on the previous arguments, we sketch three 
possible scenarios to explain the spectral state of hard 
power-law ULXs such as XI and X2, and their variability 
(or lack of) between states. 

a) Intermediate-mass BHs in the low/hard state. In 
this scenario, the hard power-law spectra are attributed 
to a true low/hard state. If low/hard to high/soft 
state transitions (that is, between a non-thermal, 
radiatively-inefficient and a thermal, radiatively-efficient 
accretion mode) must always occur for all accreting BHs 
at Lx ~ a few per cent of L^dd, this scenario implies 
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Fig. 5.— Top panel: combined 2007 XMM-Newton/EPIC spec- 
trum of XI, and \ 2 residuals. The model is an absorbed power-law. 
The best-fitting parameters are listed in Table 2. We have not in- 
cluded in the fit three "bad" channels (at pb 0.6, 0.8 and 2 keV) 
that were clearly very discrepant from any fitting model, leading to 
an additional A\ 2 ~ 30; we are unable to pin down the reason for 
the odd behavior of those three channels, but we suspect they are 
related to the EPIC-pn behavior at the edge of the chip gap (they 
are not discrepant in the MOS). Including those three channels 
does not change the best-fitting parameters. Bottom panel: the 
same spectrum cannot be well fitted by a disk-blackbody model, 
which highlights the lack of spectral curvature. 

that ULXs with a hard (T < 2) power-law spectrum 
and luminosities ~ 10 40 erg s^ 1 must be powered by 
intermediate-mass BHs with masses > a few 10 3 M Q 
(Winter, Mushotzky & Reynolds 2006). In this scenario, 
ULXs represent a different physical class of BHs; the 
reason they never switch to the high/soft state is that 
the accretion rate from the donor star in never high 
enough. However, there is still little supporting evidence 
or theoretical justification for such intermediate-mass 
BHs. Hence, we argue that this is the least likely 
interpretation, at this stage. 

b) Stellar-mass BHs in the super-Eddington state(s). 
In this scenario, power-law ULXs (with or without soft 
excess) and convex-spectra ULXs are all varieties of the 
very high state or super-Eddington regime, character- 
ized by a mix of Compton scattering, photon trapping, 
slim-disk energy advection and massive radiation-driven 
outflows (with possible transitions between substates, as 
seen for example in IC342 XI: Kubota et al. (2001)). 



This implies BH masses ~ 10Af Q , isotropic X-ray 
luminosities ~ l-30LEdd, and mass accretion rates rh ~ 
a few to ~ 1000. In this scenario, ULXs represent a 
different accretion state of normal stellar-mass BHs; 
their canonical high/soft state is not observed because 
their accretion rate is always above Eddington, and so 
is their luminosity. 

c) Massive stellar BHs without a high/soft state. In 
this scenario (somewhat intermediate between the other 
two), BHs in ULXs have masses ~ 30-IOOMq (at the 
upper end of the high-mass X-ray binary distribution), 
and isotropic luminosities Lx ~ 0.1 to a few LEdd- The 
main observational motivation for this massive-stellar- 
BH scenario comes from the observed downturn in the 
ULX luminosity distribution at Lx ~ 2 x 10 40 erg s~ 4 
(Grimm, Gilfanov & Sunyaev 2003; Swartz et al. 2004), 
which can be interpreted as ~ the characteristic Edding- 
ton luminosity for the ULX population. This implies that 
the isotropic luminosity of most ULXs is below Edding- 
ton. In turns, this means that ULXs must follow a differ- 
ent canonical sequence of state transitions than Galac- 
tic BHs, because they are always dominated by Comp- 
tonized emission at luminosities below, near and slightly 
above Eddington, without a "naked-disk" high/soft state 
in between. If there is a direct disk component, it is al- 
ways less important than the Comptonized component, 
or it may be a slim-disk component for the most luminous 
sources. Following this interpretation, there is no longer 
a gap between low/hard and very high state, just a con- 
tinuous range of coronal optical depths (getting thicker at 
higher rh) and temperatures (getting cooler); one "state" 
can morph directly into the other. 

4.3. Two- component accretion flows? 

We argue that the scenario c) above is the most con- 
sistent with multiband observations, although there may 
also be some ULXs in the highly super-Eddington state 
of scenario b). One issue that can differentiate the two 
scenarios is the relative change in rh necessary to pro- 
duce a change in luminosity. For scenario c), changes in 
luminosity between ~ 0.1 and ~ lLEdd can correspond 
to similar relative changes in the accretion rate. For sce- 
nario b), the radiative luminosity is strongly saturated 
above Eddington: L ~ (1 + mm); changes in luminos- 
ity between ~ 1 and ~ lOLEdd require an increase in 
rh by several orders of magnitude, which may severely 
constrain the mechanism responsible for the transition, 
and may have stronger detectable effects on the optical 
counterpart, or outflow density. 

Why would ULXs lack a canonical high/soft state? 
Their BHs would be a few times, or perhaps even an 
order of magnitude more massive than Galactic BHs — 
but this should not directly affect their state-transition 
properties. We suggest that the type of donor stars and 
the mechanism of mass transfer can cause a different 
spectral-state behaviour in Galactic BHs and ULXs. The 
canonical spectral state classification is based mostly on 
the appearance of Galactic BHs in low-mass X-ray bina- 
ries, also known as soft X-ray transients. Instead, the 
most luminous ULXs are thought to have OB donors, 
with masses > 10M Q , based on theoretical and observa- 
tional considerations. Stellar evolution models (Rappa- 
port, Podsiadlowski & Pfahl 2005) show that OB stars 



can provide the required mass transfer rates > 1O _6 M0 
yr _1 over their nuclear timescale ~ 10 6 yr, with a peak 
as the donor evolves to the blue-supergiant stage. And 
optical observations of ULX counterparts show that they 
are indeed predominantly located in OB associations or 
regions of recent star formation, even though it is al- 
ways difficult to identify the donor star, and separate 
the optical contribution of the donor star from that of 
the irradiated accretion disk (Ramsey et al. 2006; Grise 
et al. 2008). A typical ULX donor star is also likely 
to be filling its Roche lobe, because the fraction of gas 
captured by the BH from wind accretion alone may not 
be sufficient to produce the observed X-ray luminosity. 
We speculate that ULX donor stars may provide both 
high-angular-momentum gas, feeding the disk through 
Roche-lobe overflow, and low-angular-momentum gas, 
feeding the Comptonizing medium and jet through a sub- 
Keplerian wind, perhaps enhanced by the X-ray irradi- 
ation of the secondary. This would be analogous to the 
two-flow scenario proposed for some Galactic BHs with 
massive donors (Smith, Heindl & Swank 2002). In fact, 
it may not be a coincidence that Cyg X-l, perhaps the 
closest Galactic analog to a ULX in terms of accretion ge- 
ometry and donor star, has only rarely been seen in a true 
high/soft state (Dotani et al. 1997); instead, most of its 
state transitions are between the low/hard state (where 
it spends the majority of its time) and what should be 
properly labelled as very high state, with radio activ- 
ity and an X-ray spectrum consisting of a strong power- 
law plus thermal component (Remillard & McClintock 
2006; Feng, Tennant & Zhang 2003; Zhang et al. 1997,b). 
We note as another possible analogy that the accretion 
mode of radio galaxies (as indicated by their emission- 
line type and radio power) may depend on their source 
of gas (Bondi accretion of the hot gas and/or disk accre- 
tion of cold gas), not simply on the total accretion rate 
(Hardcastle, Evans & Croston 2007). 

The arguments above may at least qualitatively ex- 
plain why ULXs can retain a hard power-law spectrum 
even at X-ray luminosities ~ 0.1-lLEdd- The other fac- 
tor that may determine the ULX state behaviour is their 
high accretion rate, as mentioned earlier. BH transients 
can switch to a soft, disk-dominated state during an out- 
burst only if m < 1; in that case, the accretion flow 
can settle into a quasi-steady standard-disk configura- 
tion. But this cannot occur if the accretion rate dur- 
ing their outbursts stays persistently above Eddington 
(the inner part of the disk cannot survive at m > 1). 
We suggest that some transient ULXs may switch di- 
rectly from the low/hard state (at low accretion rates) 
to the slim-disk state (which produces a convex X-ray 
spectrum), or to the Compton-scattering-dominated very 
high state (which produces a cut-off power-law spec- 
trum). The standard disk may only survive or be di- 
rectly visible at large radii, and may produce the cool 
soft excess (Roberts 2007). In this scenario, the reason 
why ULXs have long-lasting phases of super-Eddington 
mass transfer rate may be due to the different mass or 
evolutionary stage of their donor stars, compared with 
the typical donor stars of Galactic BH transients. 

4.4. Thresholds for the low/hard state 

To make progress on these issues, we need to under- 
stand more about the threshold luminosity and accretion 
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Fig. 6. — Red datapoints (lower count rate): combined XMM- 
Newton/EPlC spectrum of X2 from 2007 Jun-Jul, with \ 2 fit 
residuals. Blue datapoints (higher count rate): Chandra/ ACIS 
spectrum of X2 from 2006 April 23, with x 2 fit residuals. For both 
spectra, the best-fitting model is an absorbed power- law (Table 3). 



rate at which BHs (including Galactic X-ray binaries) 
switch from the low/hard to the high/soft state, and to 
determine whether, in some cases, a system can remain 
in a hard state up to luminosities ~ 0.1-lLEdd, and then 
evolve directly towards the very high state. The scenario 
of a universal threshold between low/hard and high/soft 
state atm~ 0.01 works if the low/hard state consists of 
a radiatively- inefficient, advection-dominated flow (Esin, 
McClintock & Narayan 1997). Such flow is not expected 
to survive at higher accretion rates, when radiative cool- 
ing becomes more efficient. However, other accretion 
models may explain the persistence of a hard state up 
to higher luminosities: for example a "luminous hot ac- 
cretion flow" (Yuan & Zdziarski 2004) may be consistent 
with hard-state luminosities up to at least O.lLEdd- Fur- 
ther investigations should determine whether the hard 
state in ULXs (and perhaps in some classes of AGN) 
can extend to even higher luminosities, if the domi- 
nant power-law component is produced by a magnetized 
corona or a jet, perhaps favored by sub-Keplerian ac- 
cretion flows. Observationally, low/hard to high/soft 
transitions during several outbursts of the Galactic BHs 
GX 339-4 and XTE J1550-564 have been seen at differ- 
ent X-ray luminosities, sometimes as high as « O.lLEdd 
(Yu et al. 2007). The existence of a linear relation be- 
tween the peak luminosity in the low/hard and high/soft 
states, and the properties of the time delay between the 
two states, have been interpreted (Yu et al. 2007) as pos- 
sible evidence of the two-flow accretion geometry (sub- 
Keplerian flow responsible for the low/hard emission, and 
Keplerian flow feeding the disk). Hard-state luminosi- 
ties > 10 38 erg s _1 have also been seen in the Galac- 
tic BH candidate Cyg X-3 (Hjalmarsdotter et al. 2008), 
which has a Wolf-Rayet donor with a strong wind. Fi- 
nally, there is a group of at least 8 Galactic BHs (Brock- 
sopp, Bandyopadhyay & Fender 2004) that remain in 
the low/hard state throughout their outbursts. Under- 
standing why those BHs never switch to the high/soft 
state will also provide a useful comparison for the ULX 
behavior. 

If our suggestion that ULXs lack a canonical high/soft 
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state is correct, we should gain further insight by study- 
ing their radio behaviour. Accreting BHs have a steady 
jet (flat-spectrum radio-core emission) in the low/hard 
state; instead, they are radio-quiet in the high/soft state, 
when the jet is quenched (Fender, Belloni & Gallo 2004), 
and radio-flaring in the very high state. If some BHs 
never switch to the high/soft state, it will be interesting 
to determine whether they always have a jet, which may 
carry out a substantial fraction of the accretion power; 
perhaps the only transition is between a steady jet at 
lower luminosities and a flaring jet at higher luminosities. 
So far, radio observations have not yet been sensitive 
enough to detect point-like radio core emission in ULXs, 
and monitor their short-term variability in relation to 
their X-ray state. Only the much brighter, persistent, 



extended lobe emission (optically-thin synchrotron nebu- 
lae) has been identified in some ULXs (Miller, Mushotzky 
& Ncff 2005; Lang et al. 2007). 



We thank the referee for his/her suggestions on the 
significance of the slim-disk state, which allowed us to 
improve the Discussion. RS acknowledges support from 
a Leverhulme Early-Career Fellowship, a UK-China Fel- 
lowship for excellence, and from Tsinghua University 
(Beijing); he also thanks Prof. Shuang-Nan Zhang and 
Prof. Wcnfei Yu for their insightful suggestions and their 
hospitality at Tsinghua University and at the Shang- 
hai Observatory, respectively. GR acknowledges support 
from his NASA-G06-7102X and NNX07AR90G grants. 



REFERENCES 



Arnaud, K. A. 1996, Astronomical Data Analysis Software and 
Systems V, ASP Conference Series Vol. 101, G. H. Jacoby and 
J. Barnes eds, 17 

Begelman, M. C. 2002, ApJ, 568, L97 

Begelman, M. C. 2006, ApJ, 643, 1065 

Berghea, C. T., Weaver, K. A., Colbert, E. J. M., & Roberts, T. 

P. 2008, ApJ, 687, 471 
Brocksopp, C, Bandyopadhyay, R. M., & Fender, R. P. 2004, 

New A, 9, 249 

Chakrabarti, S. K., & Titarchuk, L. G. 1995, ApJ, 455, 623 
Chakrabarti, S. K. 1997, ApJ, 484, 313 
Done, C, & Kubota, A. 2006, MNRAS, 371, 1216 
Dotani, T., et al. 1997, ApJ, 485, L87 

Esin, A. A., McClintock, J. E., & Narayan, R. 1997, ApJ, 489, 865 
Esin, A. A., Narayan, R., Cui, W., Grove, J. E., & Zhang, S.-N. 

1998, ApJ, 505, 854 
Fender, R. P., Belloni, T. M., & Gallo, E. 2004, MNRAS, 355, 1105 
Feng, H., & Kaaret, P. 2005, ApJ, 633, 1052 

Feng, Y. X., Tennant, A. F., & Zhang, S. N. 2003, ApJ, 597, 1017 
Ferrarese, L., et al. 2000, ApJS, 128, 431 

Grimm, H.-J., Gilfanov, M., & Sunyacv, R. 2003, MNRAS, 339, 
793 

Grise, F., Pakull, M. W., Soria, R., Motch, C, Smith, I. A., Ryder, 

S. D., & Bottcher, M. 2008, A&A, 486, 151 
Haardt, F., & Maraschi, L. 1993, ApJ, 413, 507 
Hardcastle, M. J., Evans, D. A., & Croston, J. H. 2007, MNRAS, 

376, 1849 

Hjalmarsdottcr, L., Zdziarski, A. A., Larsson, S., Beckmann, V., 
McCollough, M., Hannikainen, D. C, & Vilhu, O. 2008, MNRAS, 
384, 278 

Jester, S. 2005, ApJ, 625, 667 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., 
Bajaja, E., Morras, R., & Poppel, W. G. L. 2005, A&A, 440, 775 
King, A. R. 2008, MNRAS, 385, L113 

King, A. R., Davies, M. B., Ward, M. J., Fabbiano, G., & Elvis, 

M. 2001, ApJ, 552, L109 
Komossa, S., & Schulz, H. 1998, A&A, 339, 345 
Kubota, A., Mizuno, T., Makishima, K., Fukazawa, Y., Kotoku, 

J., Ohnishi, T., & Tashiro, M. 2001, ApJ, 547, L119 
Lang, C. C, Kaaret, P., Corbel, S., & Mercer, A. 2007, ApJ, 666, 

79 

Lindblad, P. O. 1999, A&ARv, 9, 221 

Maccarone, T. J., Gallo, E., & Fender, R. 2003, MNRAS, 345, L19 
Makishima, K. 2007, in the Proceedings of the IAU Symposium 

238, Prague (Czech Republic), August 2006, V. Karas and G. 

Matt eds, Cambridge: Cambridge University Press, p. 209 
Makishima, K, Maejima, Y., Mitsuda, K., Bradt, H. V., Remillard, 

R. A., Tuohy, I. R., Hoshi, R., & Nakagawa, M. 1986, ApJ, 308, 

635 

Makishima, K., et al. 2000, ApJ, 535, 632 
Markowitz, A., & Uttley, P. 2005, ApJ, 625, L39 



Merloni, A., Heinz, S., & di Matteo, T. 2003, MNRAS, 345, 1057 
Miller, M. C, & Colbert, E. J. M. 2004, IJMPD, 13, 1 
Miller, N. A., Mushotzky, R. F., & Neff, S. G. 2005, ApJ, 623, L109 
Miyawaki, R., Makishima, K., Yamada, S., Gandhi, P., Mizuno, T., 

Kubota, A., Tsuru, T., & Matsumoto, H. 2008, PASJ, in press 
Mizuno, T., et al. 2007, PASJ, 59, 257 
Narayan, R., & Yi, I. 1994, ApJ, 428, L13 
Ohsuga, K, & Mineshige, S. 2007, ApJ, 670, 1283 
Page, M. J., Davis, S. W., & Salvi, N. J. 2003, MNRAS, 343, 1241 
Pakull, M. W., & Mirioni, L. 2002, in the unpublished proceedings 

of the symposium "New visions of the X-ray Universe in the 

XMM-Newton and Chandra era", ESTEC, The Netherlands 

(November 2001), astro-ph/0202488 
Poutanen, J., & Svensson, R. 1996, ApJ, 470, 249 
Ramsey, C. J., Williams, R. M., Gruendl, R. A., Chen, C.-H. R., 

Chu, Y.-H., & Wang, Q. D. 2006, ApJ, 641, 241 
Rappaport, S. A., Podsiadlowski, P., & Pfahl, E. 2005, MNRAS, 

356, 401 

Remillard, R. A., & McClintock, J. E. 2006, ARA&A, 44, 49 
Roberts, T. P. 2007, Ap&SS, 311, 203 
Roussel, H., et al. 2001, A&A, 369, 473 

Sandagc, A., & Tammann, G. A. 1981, Revised Shapley-Ames 
Catalog of Bright Galaxies (Carnegie Inst, of Washington) , Publ. 
635 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Smith, D. M., Heindl, W. A., & Swank, J. H. 2002, ApJ, 569, 362 

Soria, R., Baldi, A., Risaliti, G., Fabbiano, G., King, A., La Parola, 

V., & Zezas, A. 2007, MNRAS, 379, 1313 
Stobbart, A.-M., Roberts, T. P., & Wilms, J. 2006, MNRAS, 368, 

397 

Strateva, I. V., & Komossa, S. 2008, ApJ, in press 
Sunyaev, R. A., & Titarchuk, L. G. 1980, A&A, 86, 121 
Swartz, D. A., Ghosh, K. K., Tennant, A. F., & Wu, K. 2004, ApJS, 
154, 519 

Tananbaum, H., Gursky, H., Kellogg, E., Giacconi, R., & Jones, C. 

1972, ApJ, 177, L5 
Tsunoda, N., Kubota, A., Namiki, M., Sugiho, M., Kawabata, K., 

& Makishima, K. 2007, PASJ, 58, 1081 
Watarai, K., Mizuno, T., & Mineshige, S. 2001, ApJ, 549, L77 
Winter, L. M., Mushotzky, R. F., & Reynolds, C. S. 2006, 649, 730 
Yu, W., Lamb, F. K., Fender, R., & van der Klis, M. 2007, ApJ, 

663, 1309 

Yuan, F., & Zdziarski, A. A. 2004, MNRAS, 354, 953 

Zhang, S. N., Cui, W., Harmon, B. A., Paciesas, W. S., Remillard, 

R. E., & van Paradijs, J. 1997, ApJ, 477, L95 
Zhang, S. N., Mirabel, I. F., Harmon, B. A., Kroeger, R. A., 

Rodriguez, L. F., Hjellming, R. M., & Rupen, M. P. 1997, AIP 

Conference Proceedings, 410, 141 



